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ABSTRACT 

A near- wall two-equation turbulence model of the K—e type is developed for the descrip- 
tion of high-speed compressible flows. The Favre-averaged equations of motion are solved in 
conjunction with modeled transport equations for the turbulent kinetic energy and solenoidal 
dissipation wherein a variable density extension of the asymptotically consistent near-wall 
model of So and co-workers is supplemented with new dilatational models. The resulting 
compressible two-equation model is tested in the supersonic flat plate boundary layer - with 
an adiabatic wall and with wall cooling - for Mach numbers as large as 10. Direct compar- 
isons of the predictions of the new model with raw experimental data and with results from 
the K — u> model indicate that it performs well for a wide range of Mach numbers. The 
surprising finding is that the Morkovin hypothesis, where turbulent dilatational terms are 
neglected, works well at high Mach numbers provided that the near wall model is asymptoti- 
cally consistent. Instances where the model predictions deviate from the experiments appear 
to be attributable to the assumption of constant turbulent Prandtl number - a deficiency 
that will be addressed in a future paper. 
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1. INTRODUCTION 

The direct numerical simulation of compressible turbulent flows - at high Reynolds num- 
bers with all scales resolved - will not be possible for the foreseeable future, if ever at all. 
Turbulence modeling will continue to play a crucial role in the computation of high-speed 
aerodynamic flows associated with the design of advanced aircraft. In these applications, 
near-wall turbulence modeling is extremely important for the accurate prediction of wall 
transport properties, such as the skin friction and heat transfer coefficients, which are pivotal 
for design. Despite its technological importance, progress in near-wall turbulence modeling 
has been slow. Many of the commonly used near-wall models - which typically contain a 
variety of ad hoc wall damping functions - are not asymptotically consistent and yield poor 
predictions even in simple incompressible boundary layers (see Patel, Rodi and Scheuerer 1 , 
Myong and Kasagi 2 and Speziale, Abid and Anderson 3 ). These deficiencies can be fatal 
when turbulence models are applied to separated flows, or other complex boundary layers, 
that require the governing equations to be integrated directly to the wall. 

High-speed compressible flows present a whole range of new problems to near- wall turbu- 
lence modeling. Shock/boundary layer interactions with turbulence amplification and flow 
separation represent but two examples. Two-dimensional equilibrium turbulent boundary 
layers for supersonic flows are less of a problem provided that the external Mach number 
is not too large. For these flows it is generally believed that Morkovin’s hypothesis 4 - com- 
monly interpreted to mean that the turbulence statistics are only altered by compressibility 
effects through changes in the mean density - is valid for at least the range 0 < M rx> < 5. 
This hypothesis allows for the use of variable density extensions of existing incompressible 
turbulence models for which dilatational effects are neglected. However, even the ability of 
these models to reliably predict mean velocity profiles in two-dimensional equilibrium bound- 
ary layers for Mach numbers > 5 has been recently called into question (see Bradshaw, 
Launder and Lumley 5 and Huang, Bradshaw and Coakley 6 ). Of course, for non-equilibrium 
compressible flows involving shocks - or extremely high Mach numbers in or near the hyper- 
sonic flow regime - issues related to near wall turbulence modeling are even more unsettled 

(c.f. Bushnell 7 ). 

In this paper, a new near-wall two-equation turbulence model of the K — e type is 
developed for high-speed compressible flows. Two features distinguish this model from earlier 
work: (a) a variable density extension of the asymptotically consistent near-wall model of 
So et al. 8 for incompressible flows is used, and (b) high-Reynolds number models for the 
dilatational terms are implemented as developed recently by Sarkar et al. 9,10 and Coleman 
and Mansour 11 based on an analysis of homogeneous turbulence. The resulting compressible 
model is tested in supersonic flat plate boundary layers at zero pressure gradient - both 
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with an adiabatic wall and with wall cooling — for external Mach numbers in the range 
0 < Mqo < 10 . Comparisons are made with experimental data, and with the predictions of 
the K — <x) model of Wilcox 12 , in an attempt to address the following questions: 

(a) Can incompressible two-equation models be extended to supersonic boundary layers 
with a minimum of compressible corrections that are systematically determined? 

(b) Will these models perform well at high Mach numbers? 

(c) Where will the models break down, and what deficiencies give rise to it? 

These issues will be addressed in detail in the sections to follow and recommendations 
will be made for future research. 

2. THE TWO-EQUATION TURBULENCE MODEL 

We consider the supersonic turbulent flow of an ideal gas with bulk viscosity and body 
forces neglected. The equations of motion are given by: 

Mass 

^ + (P u i),i = 0 (1) 

Momentum, 

d 2 

fo(P u i) + (P u i u j),j = ~P,i ~ 3 + [p(u itj + (2) 

Energy 

d 3 

Tj^(pCpT) -f- (puiC p T) ti = -jjj + UiP'i + (TijUij + (kT t i) t i (3) 

where 

P = PRT (4) 

2 

^ij ~ ^P u k,k^ij + p{Ui t j + Ujj ) ( 5 ) 

are, respectively, the thermodynamic pressure and viscous stress tensor given that p is the 
density, u,- is the velocity, T is the temperature, R is the ideal gas constant, p is the dynamic 
viscosity, k is the thermal conductivity, and C p is the specific heat at constant pressure. Here, 
the Einstein summation convention applies to repeated indices and (•),,• denotes a gradient 
with respect to the spatial coordinate 
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Any flow variable T can be decomposed into mean and fluctuating parts in two ways: 


jr _ jc + 


( 6 ) 


where T represents the traditional ensemble average whereas T represents a Favre average 
defined by 



II 

(7) 

The mean equations of motion for compressible turbulence take the form: 



+ (wh = 0 

(8) 


0 2 

Q-pUi) + (pUiU^J = -p ti - - (pUj,j),i + [p(Ui,j + Uj,i)] tj - ( pTij ),j 

(9) 


foipCpT) + ( pUiC p T ) ti = + Uip ti + u'!p- + u\p\ { + WijUij 

(10) 


-{■(Tiju'lj + p£ [pCpQi) t i + 


where 

Tij = u'lu'l 

(11) 


Qi = uff" 

(12) 


2 

CTij = ——fAUk^Sij + H" ^j,i) 

(13) 


i> e = 

(14) 


are, respectively, the Reynolds stress tensor, the Reynolds heat flux, the mean viscous stress 
tensor, and the turbulent dissipation rate. Eqs. (8) - (10) are derived subject to one major 
assumption: turbulent fluctuations in the viscosity, thermal conductivity and specific heat 
can be neglected. 

The pressure gradient- velocity correlation can be written in the equivalent form (see 
Speziale and Sarkar 13 ) 


u 'iP',i = -{pRTu'Hi + (pRu"T")j - p'u' iti . 


(15) 


Consequently, in order to achieve closure, models are needed for: (a) the Reynolds stress ry, 
(b) the Reynolds heat flux Qi , (c) the turbulent dissipation rate e, (d) the pressure dilatation 
correlation an( l (e) the mass flux u'{. Consistent with the recent work of Sarkar et al. 9 , 
the dissipation is decomposed into solenoidal and compressible parts: 

€ = S s -\r Ec (16) 
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where for homogeneous turbulence pe 3 = and pe c = given that u)\ is the 

fluctuating vorticity. Here, e 9 represents the dissipation associated with the energy cascade. 
The length and time scales will be built up from e 3 which is to be obtained from a modeled 
transport equation. The Sarkar et al . 9 model for the compressible dissipation will be used 
in the form: 

e c = aM?e a (17) 

where M t = (u"u'-/'fRT)i is the turbulence Mach number given that 7 = C p /C v is the 
ratio of specific heats. Here, the constant a is approximately 0.5 based on direct numerical 
simulations of homogeneous shear flow 10 . 

The Reynolds stress tensor is modeled in the standard eddy viscosity form 


2 K 2 ~ 1 ~ 

Ttj = »K6 (i - 2 - ^SuSii) 

O £3 o 


( 18 ) 


where K = \th is the turbulent kinetic energy, S{j = \(dui/dxj + duj/dxi) is the Favre- 
averaged rate of strain tensor, C p is a dimensionless constant taken to be 0.096, and / M is 
a wall damping function. The wall damping function - which goes to one sufficiently 
far from the wall - will be discussed in more detail later. A simple eddy viscosity model 
was chosen since we are dealing with two-dimensional attached boundary layers; for other 
applications involving more complex flows, (18) can be easily generalized to an algebraic 
stress model or an anisotropic eddy viscosity model (see Rodi 14 and Speziale 15 ). 

The Reynolds heat flux and mass flux terms are modeled by the standard gradient trans- 
port hypotheses: 

< 19 ) 



pcr p e 3 P,z 


( 20 ) 


where Pry is the turbulent Prandtl number which is taken to be 0.9 and cr p is the mass flux 
constant which assumes a value of 0.5. A recently derived model by Sarkar 10 is used for the 
pressure dilatation correlation which is given by 


- - aipVM f + aipe 3 M? ( 21 ) 

where V = —TijUij is the turbulence production. Based on direct numerical simulations of 
homogeneous shear flow, Sarkar 10 determined that ot\ ~ 0.4 and a 2 ~ 0.2. 

In order to achieve closure, modeled transport equations for the turbulent kinetic energy 
K and the solenoidal part of the turbulent dissipation rate e 3 are needed. The exact transport 
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equation for the turbulent kinetic energy is given by 


13 


JT,{p K ) + (pUiK),i = -pTijUij - pe + p'u'i'i - u'!p- + u'/tTijj "b ^i,i "b ^i,i 


dt 


( 22 ) 


where 

V? = ~\p u i u j u j - lM, Vi 1 = u ' k a ' k i (23) 

are, respectively, the turbulent and molecular diffusion terms. If turbulent fluctuations in 
the viscosity are neglected, as well as other higher-order terms, the molecular diffusion can 
be approximated as 

T>™ = pK,i (24) 

in the boundary layer. Formally, (24) is the leading order part of (23) close to the wall 
provided that turbulent fluctuations in the viscosity and density can be neglected. The 
standard gradient transport hypothesis is applied to the modeling of the turbulent diffusion 
term: 

■Df = 12-K, j (25) 

where again px = ’pC )l f )l K 2 le a is the eddy viscosity and ctk is a constant taken to be 0.75. 
This leads to the final modeled transport equation for K : 

d 


(, pK ) + ( pUiK),i = -pTijUij - p( 1 + aM?)e s - otipVM f 


dt 


+a 2 pe s M 1 - u'!p yi + + [(m + K ,i 


(26) 


The exact transport equation for the solenoidal part of the turbulent dissipation rate is 
of the general form: 


d 

Oi(p £ °) + (pUi£s),i = 


■~pe s Ui'i + ps 3 --^ + V e , — + V>s a + (p£s,i),i (27) 


where V s , is the production of solenoidal dissipation, is the corresponding destruction 
term, and V Ss is the turbulent diffusion of solenoidal dissipation. For the sake of brevity, the 
details of the higher-order correlations that comprise V e , , and V Ss are not given (these 
correlations are quite complicated). In deriving the viscous term in (27), turbulent fluctua- 
tions in the molecular viscosity have been neglected consistent with the earlier derivations. 
The term V e , in (27) represents the production of solenoidal dissipation by deviatoric mean 
strains as well as by density fluctuations. This term will be modeled as a variable density 
extension of its incompressible form with density fluctuations neglected 13 . Hence, we take 



( 28 ) 
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where C e \ is a dimensionless constant taken to be 1.5. The same approach is used in the 
modeling of the destruction term which yields 

e 2 

$c a = C e2 f 2 pr± (29) 

where f 2 is a wall damping function and C s2 is a dimensionless constant which is taken to 
be 1.83. The standard gradient transport hypothesis is used for the turbulent diffusion term 

®.. = (^.,<) . (30) 

where cr s is a dimensionless constant which is taken to be 1.45. 

Modeling the variable viscosity term in (27) is a bit trickier. For the high Mach number 
flows to be considered in this paper, variations of the viscosity with temperature must be 
accounted for. In this vain, we use Sutherland’s law 16 for which 



where To and p Q are the reference temperature and viscosity whereas x is a constant that 
depends on the gas (for air, x = 110°/T). Eq. (31) could also be used for the viscosity trans- 
port term (1/u^Du/Dt in (27). However, this would couple the dissipation rate transport 
equation to the continuity and energy equations at the highest time derivative - an undesir- 
able feature that can cause numerical stiffness. Hence, following the work of Coleman and 
Mansour 11 , we will use the power law approximation 

M = £°(|r) , (32) 

(where n « 0.7) for the formulation of the viscosity transport term ( l/u)Du/ Dt in (27); 
Sutherland’s law will, however, be used for the calculation of the viscous diffusion terms since 
it is more accurate. Coleman and Mansour 11 showed that, for an isentropic compression, 
(32) yields the relation 

1 Du 

V~Dt = f 1 ~ n ^ ~ (33) 

where 7 = C v /C v is the ratio of specific heats. On physical grounds one would expect (33) 
to constitute the leading order part of (l/u)Du/Dt and it will be used to avoid numerical 
stiffness problems. The final form of the modeled dissipation rate equation is obtained by 
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substituting (28) - (33) into (27) which yields 



(34) 


where n « 0.7. The mean viscosity p is approximated by Sutherland’s law in (34) as well as 
in the other transport equations. 

3. NEAR WALL MODIFICATIONS FOR COMPRESSIBLE FLOWS 

Now we will address the issue of near-wall modeling in more detail. Two fundamental 
assumptions will be made: (i) the functions are analytic so that they can be expanded in a 
Taylor series near the wall, and (ii) turbulent fluctuations of the density and temperature 
vanish at the wall. While the former assumption is well accepted, the latter one - despite 
its widespread use - is somewhat debatable. However, since this assumption is rigorously 
valid for isothermal wall conditions - and a good approximation for the adiabatic wall - we 
are justified in using this simplification for the present paper. Without this assumption, the 
near wall asymptotics become much more uncertain. 

A Taylor expansion then yields: 

p' = PiV 4- fay 2 + • • • (35) 

T" = 7ij/ + 722 / 2 + • • • (36) 

where and 7 ; are functions of x,z and t given that y is the coordinate normal to the 
wall. Due to (35) and the continuity equation (1), it follows that the fluctuating velocity 
components (it, v, w) have the expansions 


u - aiy 4- a 2 y 2 + • • • 

(37) 

v — b 2 y 2 + b^y 3 + • • • 

(38) 

w = c x y + c 2 y 2 + • • • 

(39) 


near the wall based on either standard or Favre averages. Equations (37) - (39) are identical 
to their incompressible counterparts - a simplification that arises from the assumption of 
vanishing density fluctuations at the wall. By using (37) - (39), it is straightforward to show 
that we have the following near-wall asymptotics for the crucial turbulence correlations: 

K ~ 0(y 2 ), £. = 0(1), e c = 0(y 2 ) (40) 
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(41) 

(42) 

(43) 


fx T = 0(y 3 ), p'u'i i - 0(y 2 ) 

= 0(y 2 ), (<),,• = 0(y 2 ) 

V S3 =0(y), *,.=0(1) 


(the precise asymptotic behavior of the turbulent diffusion terms (25) and (30) are not 
of consequence since they constitute higher-order terms near the wall). The asymptotic 
constraints (40) - (43) can be satisfied identically with only two damping functions, and 
f 2 , which behave as 


f» = 0(y~ l ), f 2 = 0(y 2 ) 


(44) 


near the wall. 

Consistent with the underlying assumption that the density fluctuations vanish at the 
wall - which renders the fluctuating velocity to be solenoidal to the leading order near the 
wall - a variable density extension of the incompressible wall damping functions of So et 
al. 8 will be used. Validations of these wall damping functions have been carried out for a 
number of two-dimensional flows, including internal and external flows as well as flows with 
heat transfer 17 19 . In all of these cases, good correlations with direct simulation data 20-24 


and measurements 25 have been obtained. These wall damping functions are implemented in 
the form: 

U = (i + tanh(y+/115) (45) 


where 




3 U2 ( ef \ 

2 C s2 \e s ej 


(46) 


y + = yu T /i ' (47) 

f w2 = e-( R 'M\ R t = K 2 jve s (48) 

£* = £*- 2 vK/y 2 (49) 

e, = £s - 2F (VK) ti (\/K) ti (50) 


In (47) - (50), v — p/p is the local mean kinematic viscosity and u T is the friction velocity 
defined by 

u T = \/r w /p w . (51) 

where t w is the wall shear stress 26 . 
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4. THE FLAT-PLATE BOUNDARY LAYER EQUATIONS 


The steady boundary layer form of the mean turbulence equations corresponding to 
this model are now provided. Consistent with most practical computations of compressible 
flows 27 , we will solve the energy equation in its total enthalpy form where 


H = C V T + + K 


(52) 


is the Favre-averaged total enthalpy. The resulting boundary layer equations at zero pressure 
gradient take the form: 

ft r) 

(53) 


d . , d , 

+ ajW = 0 


du du d 


(/* + I‘t) 


du 

dy 


(54) 


_j>h oh a 

pu -^ +pv W = ai 


£ + *;) % + K 1 - ~k) + 0 ■ p^)] 


„du 

u ~d~y 


dK 

dy 


(55) 


dK _„dK ( 9u\ .. 3 , 

^ + = ^ \W + M? 

(— Pt\ dK 
v cr K ) dy 1 


f_du\ 

d ' 

Vdy) 

dy 


— de a _„de s 


i + n(7 - 1)] p ( 


du dv\ e a ( du s 

Tx + Yy) e ‘- C ’'K l ‘ T {Ty l 


-Cxhff + Yy 


(- . Pt\ de a 
v <r e ) dy J 


(56) 


(57) 


where p? = pC p f p K 2 /e s is the turbulent viscosity. A summary of the values of the constants 
is given as follows: C p = 0.096, C e i = 1.5, U e2 = 1.83, = 0.75, <r £ = 1.45, a p = 0.5, a = 

0.5, «i = 0.4, a 2 = 0.2, and n = 0.7. The damping functions are as specified in (45) - (46); 
consequently, the two equation model of So et al. 8 is recovered in the incompressible limit. 
The turbulent Prandtl number Pry is specified as 0.9; the molecular Prandtl number Pr is 
taken to be 0.74 for air and 0.70 for helium; and the specific heat ratio 7 is 1.4 for air and 
1.67 for helium. Sutherland’s law (31) is used to evaluate the molecular viscosity in (54) - 
(57) for air; for helium the power law quoted in Fernholz and Finley 28 is used. 
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A few comments are in order concerning the derivation of the transport equation for the 
total enthalpy H. We obtain this equation by adding the modeled transport equations for 
CpT 5 2 an d K . As an alternative approach, the exact transport equation for H can be 

modeled directly. However, we do not feel that this is a good approach since H is not a 
Galilean invariant quantity. 

Equations (53) - (57) are solved subject to the standard boundary conditions: 

u = 0, v = 0, K = 0, = 2u w (dVK/dy ) 2 (58) 

at the wall. Two different types of thermal wall boundary conditions will be considered: the 
adiabatic wall where dH/dy is zero and the cooled wall case where H is constant. At the 
edge of the boundary layer, the mean velocity and the total enthalpy are required to match 
the specified free-stream conditions. On the other hand, the turbulence quantities, K and 
e a , are assumed to be zero in the free stream. Thus formulated, the above equations with the 
appropriate boundary conditions can be solved numerically using the boundary-layer code 
developed by Anderson and Lewis 27 and modified by So et al. 8 

5. DISCUSSION OF RESULTS 

Two versions of the K — e model derived herein are used to carry out the compressible flat 
plate boundary-layer calculations to be presented in this section. One version consists of the 
full compressible model as given in Eqs. (53) - (57) and is designated as K-e model/1. The 
second version - which is designated as K — e model/2 - is a variable density extension of the 
incompressible two-equation model of Reference 8 wherein the explicit compressible terms in 
(56) and (57) are neglected along with the dK/dy term in (55). Furthermore, comparisons 
will also be made with results of the well-known K-u model of Wilcox 12 . The calculations 
of these two different versions of the new K-e model, as well as those for the K - u> 
model, will allow us to evaluate the range of validity of Morkovin’s hypothesis and assess 
the importance of having an asymptotically consistent near-wall correction for turbulence 
models. Therefore, the central questions posed in the Introduction can be addressed by 
comparing these results with well-documented experimental data spanning a wide range of 
free-stream Mach numbers and cooled wall conditions. 

The three different models mentioned above are used to calculate compressible flat plate 
boundary layers on adiabatic as well as cooled walls and the results are compared with 
measurements drawn from References 28 and 29. The calculations are carried out over the 
Mach number range 0 < < 10 for the adiabatic wall boundary condition and over the 

temperature range, 0 < Too/Tku, < 1 for the cooled wall case. Here, Moo is the free-stream 
Mach number, Too is the free-stream temperature and T aw is the adiabatic wall temperature. 
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Since only mean flow properties are available from Reference 28, comparisons are made with 
these measured quantities only and the reported skin friction coefficient, Cf = 2 t vj /~p o0 U% o 
(where Uoo is the free-stream velocity and is the free-stream mean density). In addition, 
comparisons of the calculated turbulence quantities are made with the results of the K — u> 
model. All calculations are carried out to the same momentum thickness Reynolds number 
(Re) as the measurements and the comparisons are made at these respective locations. Four 
sets of data are chosen from Reference 28. These are cases 55010504, 53011302, and 73050504 
with the adiabatic wall boundary condition and case 59020105 with a constant cooled wall 
temperature. M 0 0 for these cases are 2.244, 4.544, 10.31, and 5.29, respectively, and the 
corresponding values of Re are: 20,797; 5,532; 15,074; and 3,939. The variations of Cf with 
Moo and Too/T aw are compared with the van Driest II results given in Reference 29. Finally, 
the mean velocity profiles are compared with the van Driest law-of-the-wall for compressible 
boundary layers 30 ’ 31 and an assessment of the effects of compressibility on boundary-layer 
flows is attempted. 

The mean velocities are presented in two different forms: one in terms of wall variables, 
u + versus In y + , and another in terms of outer variables, u/Uoo versus y/S. Here, u + = 
u/u r , = yu r /V w and 8 is the boundary-layer thickness. In this way, the effects of the 
calculated u T and p on the mean velocities can be assessed. The mean temperature is 
plotted in the form of T/Too versus y/8 while all turbulence quantities, such as K and —uv, 
are normalized by u 2 r to give k + and — mJ + , respectively, for presentation. A comparison 
with the' van Driest law-of-the-wall 30 ’ 31 for compressible flows is also attempted; therefore, 
the mean velocities of two cases are plotted in terms of the compressible defined according 
to Reference 31 as = /„ + (p/p w )* du + . Since the K budgets for all cases considered are 
essentially similar, only the budget of I{ in the near-wall region for case 73050504 (M 0 0 — 
10.31 and T w /T r — 1.0) is presented. Here, the reference temperature T r is taken to be the 
recovery temperature for the cases with an adiabatic wall boundary condition and taken to 
be the adiabatic wall temperature for the cooled- wall case. 

The u + versus In y+ plots for the four cases are shown in Figures 1 and 2. In each figure, 
the calculated and measured skin friction coefficient Cf, as well as that determined from the 
van Driest II formula of Reference 29, are listed for comparison. A standard law-of-the-wall 
given by u + — (l/«)/n j/+ + B is also shown for comparison, where the von Karman constant 
k = 0.41 has been assumed. It is recognized that the intercept B should be a function of 
Moo, however, in these figures, B is taken to be 4.7. The actual value used here is not too 
important because the objective is to determine and compare the variation, if any, of the 
log-law slope with and T w jT r . The log-law is seen to hold true for all cases calculated 
using both the K — e and K —u models (Figures 1 and 2) and the value of k thus determined 
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for the three adiabatic wall cases is essentially 0.41. For the cooled- wall case, the calculated 
k is not equal to 0.41 and varies from model to model. The best agreement with data is 
given by the calculation using K — e model/2 while the smallest value of k is predicted by 
the K — u model. All models considered yield calculations of (7/ that are in agreement with 
measured data and with van Driest II values. The maximum error is less than 4% (Figures 1 
and 2). An exception is the cooled- wall case where the measurement is higher than the van 
Driest II value. According to Reference 28, the measured Cf in this case is not as accurate 
as in the other cases and this explains the discrepancy between the calculated and measured 
Cf shown in Figure 2. This comparison, therefore, shows that the additional compressible 
terms in the governing equations do not significantly affect the calculated results so long 
as the near-wall model used is the same and is asymptotically consistent. In other words, 
Morkovin’s hypothesis is valid for flows with free-stream Mach numbers as high as 10 and 
wall temperature ratios noticeably smaller than one. 

In the past, velocity profiles in wall coordinates were invariably plotted in terms of 
to illustrate the existence of the van Driest log-law and the constancy of k in compressible 
boundary-layer flows 6 ’ 30 ’ 31 . Since then, the compressible law of the wall is typically taken 
to be given by rather than by « + , and k is considered to be about 0.41 and constant 
over the Mach number range of 0 < < 5. The calculated and measured velocity plots 

given in Figures 1 and 2 show support for the compressible law-of-the-wall when it is written 
in terms of u + rather than uf. Furthermore, k is determined to be approximately 0.41 
and is relatively constant over the Mach number range of 0 - 10 for the adiabatic wall 
boundary condition. These results seem to conflict with the proposal of van Driest 30 . In 
order to resolve this seeming contradiction, the velocity plots of uf versus In for cases 
55010504 (Moo = 2.244, T w /T r = 1) and 53011302 (M^ = 4.544, T w /T r = 1) are shown in 
Figure 3. In addition, the compressible law-of-the-wall as given in Reference 31 is shown for 
comparison. It can be seen that a line that is parallel to the compressible law-of-the-wall 
can be drawn through a few of the data points spanning over a narrow range of j/+. On the 
other hand, the calculated profiles are in agreement with data over a wider range of y+. The 
slopes of the calculated profiles are roughly parallel to that determined from measurements 
and are slightly larger than the slope of the compressible law-of-the-wall shown. Therefore, 
irrespective of how the velocity profiles are plotted, the calculations are in good agreement 
with data. However, the slope of the log-law appears to be given by (0.41) -1 only when the 
profiles are plotted in terms of u + . 

Based on the above comparisons, it seems that there is very little difference between 
the predictions of the K — e and K — u> models. This is particularly true in regard to the 
calculations of Cf. The semi-log plots shown in Figure 1 tend to mask the differences found 
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between the models in the calculations of the mean temperature and density profiles. These 
differences begin to show up in the plots shown in Figure 3 and in the predictions of the 
cooled- wall case (Figure 2). In order to determine the actual difference between the mean 
flow predictions of the K—e and K —ui models, plots of u/U^ versus y/6 are given in Figures 
4 and 5 for the adiabatic and cooled wall conditions, respectively, while the corresponding 
mean temperature profiles are plotted in Figures 6 and 7. These results show that there 
are substantial discrepancies between measurements and the profiles of u/Uoo and T /T^ 
calculated using the K — u model. The discrepancies increase as increases and as T w /T r 
decreases. The best results are given by K — e model/1. On the other hand, the predictions 
of K — e model/2 are very close to those of K — e model/ 1 and are substantially different 
from those of the K — u model. This is further evidence that an asymptotically consistent 
near-wall model is more important than the additional compressible terms as far as the 
predictions of the mean flow properties are concerned. 

The ability of the K — u model to predict the variations of (7/ with M <*, and T w /T aw is 
well established 12 . If the proposed K — e model is to be accepted, its ability to predict (7/ 
for different M and T w /T aw has to be demonstrated. In Figure 8, the variation of <7//((7/) t - 
with Moo for the case of the adiabatic wall boundary condition is shown. Here, ( C/)i is the 
skin friction coefficient for an incompressible flow evaluated at Rg — 10 4 which is determined 
to be 2.70 x 10~ 3 . This figure shows a comparison of the calculations of K — e model/1 and 
K — e model/2 with the van Driest II curve 29 . Essentially, there is no difference between the 
predictions of K — e model/ 1 and K — e model/2 and both results are in excellent agreement 
with data. The predictions for the cooled wall case are shown in Figure 9. Calculations for 
this case are carried out for M<x> = 5 and Rg — 10 4 where the incompressible Cj is again 
determined to be 2.70 x 10~ 3 . Three sets of calculations are presented. These are for K — e 
model/1, K—e model/2 and a third version of K—e model/1 where the dK/dy term in (56) is 
neglected. It can be seen that an error of 5% or larger starts to accumulate at approximately 
T w /T aw = 0.4 for K — e. model/1. This trend is contrary to the results reported in Reference 
12. An examination of the governing equations solved by other researchers revealed that, 
besides the differences noted in the turbulence model equations, the mean enthalpy equation 
solved by these researchers does not include the term dK/dy on the right hand side of (55). 
Indeed, when the dK/dy term is neglected, an overall important improvement is obtained. 
The predicted (7/ at T w /T aw = 0.2 is increased by about 6%, thus giving a much better 
agreement with the van Driest II curve 29 . If the additional compressible terms in the K — e 
equations are further neglected (K - e model/2), the calculated (7/ is only improved by 
about 2%. The remaining disagreement could be attributed to the assumption of a constant 
turbulent Prandtl number which cannot properly account for a reduction of turbulent mixing 
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resulting from a cooled wall (in fact, the dK/dy term can be formally dropped if Pr t — 0.74; 
such a lowering of the turbulent Prandtl number would not be inconsistent with experiments 
for the cooled wall case). However, this effect can be appropriately incorporated in a heat- 
flux model and modifications can be formulated to describe its near- wall behavior 19 ’ 32 . In 
other words, if highly cooled-wall flows are to be predicted correctly, turbulent heat fluxes 
and their near-wall behavior need to be modeled. 

The distributions of k + and —uv 1 across the boundary layers are compared in Figures 
10 - 13. In these figures, only the results of I< — e model/1 and the I< — u> model are 
compared. The results for the adiabatic wall boundary condition are shown in Figures 10 
and 12 while those for the cooled-wall case are plotted in Figures 11 and 13. It can be seen 
that the K — u> model underpredicts k + in the near- wall region and overestimates its value 
in the outer part of the boundary layer. This incorrect prediction is common for all cases 
studied. The predicted near-wall behavior is substantially different from that given by K — e 
model/1. Instead of yielding a finite slope for k + at the wall, a near zero slope is calculated. 
Furthermore, the maximum k + calculated is about half that given by K-e model/ 1. On the 
other hand, the K — u> model yields the correct near- wall behavior for —uv + but overpredicts 
its value in the outer part of the boundary layer. The overprediction extends across the range, 
0.2 < y/S < 1.0. It is now clear that the K — ui model is formulated to give correct results for 
the mean velocity and the wall shear stress; however, its predictions of other properties are 
in doubt, particularly, in regard to the near-wall behavior. Reduction of turbulence activity 
in the outer part of the boundary layer is clearly evident when either compressibility or wall 
cooling effects are present. The reduction increases as M increases and T w /T r decreases. 
This is further substantiated by the very significant drop in the maximum value of k + as 
increases. Therefore, it is expected that turbulence activity will be substantially reduced in 
a flow where the free-stream Mach number is large and the wall is highly cooled. 

The predictions of the near-wall flow can be further examined. With the assumption 
of vanishing p' and T" at the wall, Taylor expansions of the fluctuating velocities, density 
and temperature about y = 0 are given by (35) - (39). These expansions together with 
the definition for e s can be rearranged to give the following dimensionless expansions for 
k + ,— uv+ and e+, i.e., 

k + = a k (y+) 2 + b k (y+) 3 + ..., ( 59 ) 

-uv* = a uv (y+f + b uv (y+) 4 + , . . , (60) 

e+ = 2a* + 4%+ + ..., (61) 

where the a’s and 6’s are time-averaged coefficients that are functions of x and 2 . From 
these expansions, it can be easily deduced that h + /ef(y+) 2 = 0.5 at the wall. Therefore, the 
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asymptotic behavior of k + /ef{y+) 2 is 0.5 and is independent of Moo as well as the thermal 
wall boundary conditions. The accuracy to which a model can predict this quantity is a 
reflection of its asymptotic consistency. Finally, the variation of ef at the wall with Mach 
number can be deduced from a plot of ak versus Moo ; therefore, this result is shown in 
Figure 14 for the cases with the adiabatic wall boundary condition. The value of ak for the 
corresponding incompressible flow is taken from Reference 8. Model calculations of ak,a uv 
and k + /ef(y+) 2 are tabulated in Table 1 for comparison. Since the K — u> model is not an 
asymptotically consistent near-wall model, its predictions of these limiting values are poor; 
therefore, they are not listed in Table 1. From these results, it can be seen that the present 
K — e model is asymptotically consistent up to a Mach number of 10 for both adiabatic 
and cooled wall conditions. In other words, the predictions of k + and — uv 1 shown in Figures 
10 - 13 are more likely to be correct compared to those given by the K — oj model. These 
results further show that a k is a function of the Mach number and decreases with increasing 
(see Figure 14). This means that viscous dissipation at the wall also decreases with 

increasing Moo. 

The near-wall budget of K for the = 10.31 case (73050504) is shown in Figure 
15. It is evident that the budget of K bears a strong resemblence to that calculated for 
incompressible flows'^. The additional compressible terms have a negligible effect on the 
near-wall K budget. Therefore, the assumptions made to extend the near-wall damping 
function f 2 in the dissipation-rate equation to compressible flows appear to be justified. 
Again, it can be seen that viscous diffusion balances dissipation at the wall. This balance 
extends to about = 4 where turbulent diffusion and production become important. In the 
region, 4 < < 15, viscous and turbulent diffusion as well as production and dissipation 

are equally important. Beyond y+ = 15, production and dissipation are approximately in 
balance, just as in the case of incompressible flows. Consequently, the near-wall behavior of K 
is very similar for both incompressible and compressible flows (thus, explicit compressibility 
effects are not that important in the near- wall region). 

6. CONCLUSIONS 

In contrast to the conclusions drawn by Bradshaw et al. 5 and Huang et al. 6 , the present 
investigation shows that conventional K — e models can be extended to calculate equilib- 
rium compressible boundary-layers, at high Mach numbers, if the near-wall modifications to 
these models are asymptotically correct and internally consistent. When properly modified, 
K - e type models can be used to calculate equilibrium compressible boundary-layers with 
free-stream Mach numbers as high as 10 and wall temperature ratios as low as 0.2. Two 
different versions of the K - e model proposed in this study are tested: one with all of 
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the systematically derived dilatational terms included in the governing equations (I< - e 
model/1) and another which is a variable density extension of the incompressible limit of 
this model where explicit compressible terms are neglected (K - e model/2). Calculated 
mean flow properties and wall shear stresses from these two versions of the I< - e model are 
in excellent agreement with measurements. The near- wall behavior of the calculated tur- 
bulence properties is consistent with exact results and the correct limiting values for these 
properties are recovered. Compressibility effects are found to have a negative influence on 
turbulent mixing which is reflected in a reduction of the maximum value of the turbulent 
kinetic energy and the value of the viscous dissipation at the wall as the free-stream Mach 
number increases. Wall cooling also produces the same effects and the two versions of the 
proposed K - e model mimic the trend very well. This means that Morkovin’s hypothesis 
is valid for the full range of Mach numbers and wall temperatures tested, provided that the 
near- wall two-equation models used to close the governing equations are asymptotically cor- 
rect and internally consistent in the near-wall region. On the other hand, the well accepted 
K ~ oj model provides good correlations with measurements for the mean velocity and wall 
shear stress only. Its predictions of the temperature profiles are in considerable error, partic- 
ularly for high free-stream Mach numbers and for the cooled-wall boundary condition. The 
reason for these discrepancies is traced to the model’s inability to reproduce the near-wall 
turbulence properly; the modeled asymptotic behavior of the turbulence properties is not 
consistent with the exact equations governing the transport of these quantities. Despite 
these deficiencies, it is still rather surprising how well the K — u> model predicts the mean 
velocity and skin friction at high Mach numbers. 

Considering the excellent results obtained in this paper for equilibrium turbulent bound- 
ary layers, future applications of this compressible K — e model are planned for the study of 
more complex non-equilibrium boundary-layer flows involving shocks. In these applications, 
we would expect the explicit compressible terms appearing in the modeled transport equa- 
tions for K and e 3 to play a more important role. After these further tests are completed, 
we will have a much better idea of the full range of applicability of this new compressible 
K — e model for wall-bounded turbulent flows. 
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Table X. Asymptotic near- wall behavior of the calculated turbulence properties. 
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Figure 1. Semi-log plots of u + for adiabatic wall boundary condition: (a) M, 
Moo = 4.544, and (c) = 2.244. 
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igure 4. Linear plots of u/U 0 0 for adiabatic wall boundary condition: (a) 
Moo = 4.544, and (c) = 2.244. 
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Figure 5. Linear plots of u/U 0 0 for cooled wall boundary condition. 
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Figure 8. Variation of Cf/(C f )i with for adiabatic wall boundary condition. 
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Figure 9. Variation of Cj/(C/)i with T w /T aw for = 5.0. 
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Figure 10. Behavior of k + across the boundary layer for adiabatic wall boundary condition: 
(a) Moo = 10.31, (b) Moo = 4.544, and (c) M^ = 2.244. 







Figure 12. Behavior of — uv 1 across the boundary layer for adiabatic wall boundary condi 
tion: (a) M 0 0 = 10.31, (b) M, ^ = 4.544, and (c) M 0 0 = 2.244. 
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